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Abstract. - The non-equilibrium ageing behaviour of the three-dimensional and four- dimensional 
critical Ising spin glass is studied for both binary and gaussian disorder. The same phenomenol- 
ogy of the time-dependent scaling as in non-disordered magnets is found but the non-equilibrium 
exponents and the universal limit fluctuation-dissipation ratio depend on the distribution of 
the coupling constants. 



Understanding the kinetics of spin systems or glasses after a rapid quench from an initial 
disordered state continues to pose challenging problems. A key insight has been the obser- 
vation that many of the apparently erratic and history-dependent properties of such systems 
can be organized in terms of a simple scaling picture [1]. The simplest way to interpret this 
is by assuming the existence of a single time-dependent length-scale L(t) in the problem. 

For simple magnets without disorder, one typically finds a power-law scaling L(t) ~ t 1 ^, 
where z is the dynamical exponent. It has been realized that the ageing behaviour is more fully 
revealed in observables such as the two-time autocorrelation function C(t, s) := {<p(t)<p(s)) or 
the two-time linear autoresponse function R(t,s) :— S((f>(t)) /Sh(s)\ h=0 , where <p(t) denotes 
the time-dependent order-parameter, h(s) is the time-dependent conjugate magnetic field, t 
is referred to as observation time and s as waiting time. One says that the system undergoes 
ageing if C or R depend on both t and s and not merely on the difference r = t — s. For 
simple magnets, these two-time functions are expected to show dynamical scaling in the ageing 
regime t,s ^>> t m i C ro and t — s ^> £ m i cr0 j where i m i cr o is some microscopic time scale. Then 



where a, b are non-equilibrium exponents. The scaling functions fc,n(y) should satisfy the 
following asymptotic behaviour 



C(t,s) =s- b f c {t/s) , R{t,s) 
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as y — y oo and where Ac and A/j, respectively, are known as the autocorrelation [2,3] and 
autoresponse exponents [4]. In what follows, we shall concentrate on quenches onto the critical 
point. Then for simple ferromagnets a = b = 2(5/ vz — (d — 2 + rj)/z where /3, rj, v are standard 
equilibrium critical exponents. Furthermore, for spatially short-ranged initial correlations, one 
has Ac = Xr. These expectations on critical ageing have been confirmed in a large variety 
of models through simulations, exact solution or field-theory calculations, see [5-8] for recent 
reviews. Critical ageing JIJ also occurs in systems without detailed balance such as the contact 
process, but now with the exponent relations 1 + a = b = 2(5 /v^z while the autocorrelation 
and autoresponse exponents still coincide [9,10]. 

It is natural to try and see whether these results might be extended to glassy systems. 
In this letter, we shall concentrate on the Ising spin glass, with a static Hamiltonian TL = 
~ j) Ji,j°'i' J j- Here <Xj = ±1 are the usual Ising spins and the nearest-neighbour couplings 
Jij are random variables. We shall consider (i) a binary distribution where Jjj = ±1 and 
(ii) a gaussian distribution of the Ji^ with zero mean and variance one. The dynamics of 
the model is given by a master equation where the rates are chosen according to heat-bath 
dynamics. It is now established [11] that this model undergoes in d > 2 dimensions an 
equilibrium phase-transition between a paramagnetic and a frustrated spin-glass phase. On 
the other hand, when considering a quench below the spin-glass critical temperature T c , there 
has been considerable debate on the precise relationship between the relevant time and length 
scales. It has been attempted to summarize the present state of knowledge into the form [12] 

i(i) ^« p ^(_A_)") (3) 

where Ao is an energy scale of order T c , if) is a barrier exponent and £(T) is the equilibrium 
correlation length at temperature T. This form has been used to fit successfuly simulational 
data in the 3D and AD Edwards- Anderson model [12,13]. Although the typical length scales 
are merely of the order of a few lattice sizes, see e.g. [14], the relaxation times are sufficiently 
large for a dynamical scaling to set in. We shall return to this below. While this expression, 
if correct, points towards a cross-over behaviour between a simple power-law scaling and an 
exponential scaling as would follow from the droplet model, it also suggests that at criticality, 
simple power-law scaling should prevail. The determination of the exponent z in spin glasses 
and in real materials is a much-studied topic, see [11] and references therein. It is the aim of 
this letter to further test this idea in the critical Ising spin-glass, paying special attention to 
the dynamical scaling behaviour of two-time functions in the ageing regime. A similar analysis 
was recently performed in a damage-spreading context [15]. From now on we consider the 
three- and the four-dimensional Ising spin glass quenched to T = T c (the assumed values 
of T c are listed in table [IJ from a fully disordered initial state. We shall study the scaling 
behaviour of the magnetic autocorrelation C(t, s) = N^ 1 J2i { (J i{t)' J i{ s )) an d of its associated 
linear response with respect to an external magnetic field h. As usual, since response functions 
are too noisy to be measured directly, we study instead the thermoremanent magnetization 
by turning on the magnetic field immediately after the quench at time t = and keeping it 
until after the waiting time s has elapsed. The magnetization measured at a later time t is 
related to the linear response function by M(t,s) = h J^du R(t,u). The systems simulated 
contained 50 3 and 20 4 spins, respectively. Some other system sizes were also briefly considered 
in order to check against finite-size effects. The autocorrelation data discussed in the following 
have been obtained after averaging over typically a few thousand different bond distributions. 
For the thermoremanent magnetization we averaged over at least ten thousand samples with 
different realizations of the couplings. 
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Fig. 1 - Scaling of (a) the autocorrelation function C(t, s) and (b) the thermoremanent magnetization 
M[t, s) of the 3D binary Edwards- Anderson spin glass at criticality. The full curves are C(t, s) = cos~ b 
and M(t, s) = mos~ a , where co, mo were fitted to the numerical data. Here and in the following error 
bars are smaller than the symbols sizes, (c) Temporal evolution of the autocorrelation C(t, 0) in the 
three-dimensional cases (b: binary, g: gaussian). 



We first study the scaling as a function of the waiting time s, with t/s fixed. In figure^we 
show data for the 3D binary case. We observe that both correlation (figure^) and integrated 
response (figure \T]p) are consistent with a power-law scaling. Similar results were obtained 
also in AD and for the gaussian case. In table [I] our results for the exponents a and b are 
listed. Their values are close to results inside the spin-glass phase but near to T c [16]. 

We observe that a = b within our numerical errors for a fixed choice of the distribution of 
the couplings. However, the results for binary and gaussian distributions are different. 

In order to check that the results of table U relate to the true long-time scaling regime, 
we compare with the expected relation a = b = l3 EA /{vz) = (d — 2 + t/ba)/(2z) where 
/3ea,t?ea are the equilibrium critical exponents of the Edwards- Anderson order parameter 
which is quadratic in the magnetization. Their values were measured many times. For the 
3D binary case we quote 7 lEA = -0.225(25) and z = 5.65(15) [17] or ?7ea = -0.337(15) [18], 
leading to a = 0.068(4) and 0.059(3), respectively. For the 3D gaussian case tjea = —0.42(3) 
and z = 6.45(10) [19,20] leading to a = 0.044(3) whereas the result r/ EA = -0.36(6) [21] 
gives a = 0.049(5). In four dimensions one has for the binary case ijea = —0.31(1) and 
z = 4.45(10) [22] leading to a = 0.19(1). Finally, for the AD gaussian case r] EA = -0.35(5) [23] 
or t^ea = -0.44(2) [19,20] and z = 4.9(4) [19,20] which gives a = 0.17(1) or a = 0.16(1). 
These results for a agree very well with the values extracted from our dynamical simulations 
in three dimensions and for the AD binary case. For the AD gaussian case our value for a is 
substantially larger than the value obtained when combining the literature values of 77EA and 
z. The origin of this discrepancy in not clear to us. 

A second non-equilibrium critical exponent is given by the power-law decay of the auto- 
correlation for long times. As shown in figure^ the quantity C(t, 0) displays a power-law 
behaviour over more than three decades, making a very precise determination of the exponent 
\cl z possible, see [3], but which does depend on the distribution of the J^-. We have checked 
that the distribution-dependence of Ac does not simply result from inaccurate determinations 
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Table I - Critical temperature and universal nonequilibrium quantities of the critical Ising spin glass, 
for both binary and gaussian distributions of the nearest-neighbour couplings. In four dimensions even 
our longest runs did not permit us to reliably determine Xr/z. 





binary 


d 


T c 


a 


b Xc/z 


Xr/z 




3 


1.19 


[17] 


0.060(4) 


0.056(3) 0.362(5) 


0.38(2) 


0.12(1) 


4 


2.0 


[24] 


0.18(1) 


0.17(1) 0.615(10) 




0.19(1) 




gaussian 


d 


T c 


a 


b Xc/z 


Xr/z 




3 


0.92 


[19] 


0.044(1) 


0.043(1) 0.320(5) 


0.33(2) 


0.09(1) 


4 


1.78 


[19] 


0.22(1) 


0.23(1) 0.68(1) 




0.16(1) 



of T c . In 3D, using T c = 1.14 [18] we find Xc/z — 0.348(5) for the binary case and using 
T c = 0.95 [21], we find Xc/z = 0.335(5) for the gaussian case. These values still appear to be 
significantly different from each other. 

Summarizing, we have found evidence that the dynamic universality class may depend on 
the choice of the distribution of the coupling constants. For the equilibrium transition, even 
the very existence of a phase-transition in 3D was questioned for a long time, see [11] and the 
question of its universality remains controversial. For example, numerical and experimental 
evidence against universality was presented in [19, 20] while a recent high-temperature study 
of the Edwards- Anderson susceptibility with four symmetric random distributions asserted 
the universality of the exponent 7 [25] in four to eight dimensions. 

Next, we turn to the form of the scaling functions themselves. In figure |3 we show the 
two-time scaling of the spin-spin autocorrelator, for several waiting times s. It is clear that 
C(t,s) does not merely depend on the time difference t — s and therefore the system ages. 
We find in all cases a nice collapse of the rescaled autocorrelator s b C(t, s) compatible with a 
simple power-law scaling L(t) ~ t x / z . 

In a similar way, in figure [3] we show the two-time scaling of the integrated response. 
Again, we observe a very nice collapse of the data in terms of a simple power-law scaling. The 




Fig. 2 - Scaling of the autocorrelation in the critical Edwards- Anderson spin-glass with the following 
dimensions and couplings: (a) 3D binary, (b) 3D gaussian, (c) 4D binary and (d) 4D gaussian. 
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Fig. 3 - Scaling of the thermoremanent magnetization in the Edwards- Anderson spin glass at criticality 
with binary couplings in (a) 3D and (b) 4D. The full curve is the prediction Jj>) of local scale- 
invariance, see text. The results of very longs runs for a single waiting time are shown in (c) 3D with 
s = 100 and (d) 4D with s = 25. 



fact that both the autocorrelator and the thermoremanent magnetization can be described 
in terms of such a power-law scaling is evidence in favour of the time-dependent length-scale 
JsJJ. Closer inspection of figure|2td, however, reveals an unexpected subtlety. In principle, one 
would like to extract an exponent X' R /z from the slopes in that figure. It turns out that the 
values of X' R /z thus obtained are significantly different (3D binary: 0.45, 3D gaussian: 0.41, 
4D binary: 0.72, 4D gaussian: 0.76) from the ones found before for \c/z. Indeed, if one goes 
to larger values of y — t/s as is shown in figure|2td, we find that our data are systematically 
above the asymptotic power-law ~ (t/ s)~ Xr I z which we obtained from smaller values of t/s. 
This passage from an effective exponent X' R / z at intermediate values of the scaling variable y 
to the truly asymptotic value Xr/z at larger values of y has also been observed in the critical 
ageing of the ferromagnetic Ising model with Kawasaki dynamcis [26,27]. We shall return to 
a more quantitative discussion of M(t, s) below. 

In simple ferromagnets, ageing occurs out of equilibrium. The distance from the equilib- 
rium state may be measured through the fluctuation-dissipation ratio 

*<<■■'> = ras;=*< f/ »> (4) 

In particular, the limit fluctuation-dissipation ratio Xoo = lim^^oo X(x) is expected to be 
an universal number [6] . In figure 01 we display the fluctuation-dissipation relation between 
the scaled autocorrelation and thermoremanent magnetization in both 3D and 4D. For small 
time differences t — s, the autocorrelation C(t,s) is large and we are in a quasiequilibrium 
regime with a fluctuation-dissipation ratio X(t,s) w 1. On the other hand, for well-separated 
times s and t, the systems moves out of equilibrium and we can read off the limit fluctuation- 
dissipation ratio X x . We remark that since we had to use the scaled forms s a C(t,s) and 
s a M(t, s) rather than the unsealed observables, X(t, s) cannot merely depend on the value of 
the autocorrelation, in distinction with what occurs in the mean-field theory of spin glasses. 
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Fig. 4 - Determination of the fluctuation-dissipation ratio in the critical Edwards- Anderson spin glass 
with binary couplings in (a) 3D and (b) 4D. The dashed lines have slope one and show the quasi- 
equilibrium behaviour while the slopes of the full lines give the limit fluctuation-dissipation ratio Xoo . 
The insets show the function X := T c M(t, s)/(hC(t, s)) as a function of s/t for both the binary and 
the gaussian cases. In (c) we show X in 3D for the choices of the distribution and of T c : binary with 
T c = 1.19 (black) and T c = 1.14 (gray), gaussian with T c = 0.92 (black) and T c = 0.95 (gray). 

Finally, in figure^ we display the ratio T c M(t, s)/(hC(t, s)) in 3D as a function of s/t and 
for several choices of T c . For widely separated times, s/t — » and this ratio should converge to 
the limit fluctuation-dissipation ratio X^. We therefore find (see tableQ} that apparently 
depends on the distribution of the couplings as well, and independently of the assumed T c . 

The qualitative dynamical scaling behaviour of the magnetic two-time observables is quite 
analogous to the one of the conceptually simpler magnets, up to modified values of the crit- 
ical nonequilibrium exponents. We now inquire into the form of the scaling function of the 
thcrmoremanent magnetization. It is known that by extending dynamical scaling eq. to a 
local scale-invariance with infinitesimal local scale-transformations 1 1— ► (1 +e) z t, r i— ► (1 +e)r 
(with an infinitesimal e = e{t, r) which depends on both time and space), the form of f R can 
be found from the requirement of covariance of i?(i, s) under these transformations, leading 
to [28,29] /r(j/) ~ y 1 + a ^ A «/ z (y — i) a . This prediction has been reproduced in many spin 
systems quenched to a temperature T < T c and whose dynamics is described by a master 
equation [7,10,28-32]. Hence, we expect M(t,s) = M s~ a fM(t/s), where 

/m(w) = V^ njz 2F 1 (1 + a,-a + \ R / z; 1 - a + \ R / z; l/y) (5) 

with a normalization constant Mq. At T = T c , however, a second-order e-expansion produces a 
deviation from that prediction [8]. As discussed previously we observe in critical spin glasses 
a passage from an effective exponents A^/z at intermediate values of y = t/s to the truly 
asymptotic value for larger y. Hence we cannot expect equation (0 to describe the scaling 
function Jm{v) for all values of y = t/s. In figure|2tb we compare the numerical data, in both 
3D and 4D with the prediction JSJ where we have inserted the values of the exponents a and 
\' R / z which have been determined earlier. While we find a nice agreement of the prediction 
of local scale-invariance with our data for y — t/s not too large, for very large arguments 
the precise behaviour of the scaling function /j\/(y) cannot be fully reproduced. 
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In summary, we have studied the dynamical scaling behaviour of critical three- and four- 
dimensional Ising spin glasses in the ageing regime. We find for these disordered systems 
evidence of a dynamical power-law scaling analogous to non-disordered critical magnets. The 
measured values of several independent universal non-equilibrium critical quantities suggest 
that the dynamic universality class may depend on the choice of the distribution of the random 
coupling constants. 
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